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When  it  came  to  selecting  a  thesis  topic,  my  hope  was  to 
do  something  related  to  or  about  the  Space  Shuttle.  I  was 
fortunate  that  in  January  1982,  Capt  Jim  Hodge  came  to  teach 
at  AFIT  from  the  Air  Force  Flight  Test  Center  with  many 
unanswered  questions  about  the  space  Shuttle  Orbiter  reentry 
heating  affects  which  may  be  of  great  concern  to  future  Air 
Force  and  NASA  mission  planning.  Despite  my  lack  of  knowledge 
in  several  areas  concerning  this  project,  instruction  and 
assistance  from  Capt  Hodge  on  the  use  of  the  AFFTC  HtA'i't.ST 
program,  from  Dr.  Will  Hankey  in  the  fundamentals  ot  high 
speed  aerodynamics,  and  from  Capt  David  Audley  in  th«? 
fundamentals  ot  Systems  Identification  theory  made  completion 
ot  this  thesis  possible.  My  special  thanks  to  them,  to  t  ho 
rest  of  the  Engineering  and  Math  department  faculty  at  amt, 
to  my  classmates,  to  my  friends  and  family  tor  their  spiritual 
support,  to  the  Air  Force  for  permitting  me  to  be  here,  and  to 
God  who  receives  all  the  glory  tor  the  things  He  has  done  for 
me . 

I  sincerely  hope  this  will  be  ot  benefit  to  lutui.r 
investigators  ot  the  Orbiter  reentry  aerodynamics. 


John  K.  Wood 
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The  Space  Shuttle  Orbiter  has  conducted  various  flight 
test  maneuvers  during  the  Orbital  Flight  Test  Program  and 
beyond  in  order  to  establish  heating  rates  at  various  angles 
ot  attack  and  center  of  gravity  positions.  The  objective  or 
this  project  was  to  investigate  the  heating  effects  on  the 
body  flap  during  the  entry  flight  test  maneuvers  and 
determine  the  limits  on  the  flap  deflections  tor  an  aft  center 
ot  gravity  envelope  on  the  Orbiter. 

An  analysis  program  available  from  the  Air  Force  flight. 
Test  Center  (AFFTC)  uses  applicable  simulation  equations  tor 
the  aerodynamic  performance  of  the  TPS  in  their  flight 
simulator  and  data  reduction  program  (HEATEST).  The  approach 
was  to  use  the  HEATEST  program  to  determine  the  heating  on  one 
critical  point  on  the  body  flap  during  the  STS-2  entry  flight 
test  inputs  (Flap  Thermocouple  IV07T95Q8).  From  these  results 
and  other  results  on  the  outboard  elevons,  the  maximum  flap 
deflection  angle  and  Orbiter  aft  center  of  gravity  eoulu  tit- 
established. 

Using  the  HEATEST  program,  results  were  obtained  Lot 
angle  ot  attack  and  body  flap  deflection  effects  but 
difficulty  was  encountered  with  Reynold's  number  el  feet. 
Several  theoretical  approaches  were  investigated  resulting  in 
at  least  one  technique  that  could  be  applied  to  HEATEST  i-; 
linearize  or  eliminate  the  Reynold's  numb,  r  effect. 
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L. _ Introduction 


Background 

The  United  States  Space  Transportation  system,  the  Space 
Shuttle,  has  presented  ae rody nam i c i s t s  with  a  unique  new 
challenge.  Development  ot  the  entry  flight  cha  rac  te  r  i  s  t:  i  cs 
and  range  ot  operations  ot  its  primary  component,  the  orbiter, 
seen  returning  in  Fig.  1,  is  vital  to  the  overall  success  of 
its  mission.  A  systematic  approach  ot  integrating  the  best  or 
ground  tests,  flight  simulat ions,  and  flight  tests  needed  t<> 
be  developed. 

Xtlfl-  AFFTC  Role 

The  Air  Force  Flight  Test  Center  (AFFTC)  has  been  and  is 
continuing  to  conduct  an  evaluation  or  the  Orbiter  to  assess 
its  capability  to  meet  Department  of  Defense  requirements. 
The  goal  ot  the  evaluation  is  to  identity  and  remove  placards 
in  the  Orbiter's  operational  envelope  that  could  uttect 
nominal  and  contingency  landing  operations  trom  vanoenbery  Air 
Force  Base,  California,  which  are  more  severe  than  the  p re.. eat 
launch  and  landing  mode. 

AFFTC  Models 

Equations  and  parameters  were  chosen  to 
aerodynamic  heating  and  one-di  mens  iona  1  thermal  models  i  i  . 
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AFFTC  flight  simulator  (Ref  1).  The  simulator  was  also  used 
for  flight  planning  and  the  design  of  flight  test  maneuvers  to 
enhance  flight  test  data  reduction  tor  angle  of  attack  and 
center  of  gravity  envelope  expansion.  Rapid  expansion  of  the 
flight  envelope  is  required  prior  to  more  severe  entry 
conditions  and  prior  to  removal  of  the  flight  test  data 
instrumentation  from  the  operational  vehicles.  Upon  obtaining 
embedded  thermocouple  data,  a  reduction  technique  based  on 
systems  identification  theory  was  employed  (Ret  A). 

In  systems  identification  theory,  models  are  simulation 
equations  and  parameters  required  tor  aerothermody num 1 o 
performance  of  the  arbiter's  Thermal  Protection  System  (TPs). 
The  AF FTC  simulator  data  base  goes  beyond  the  simple  use  or  a 
sphere  (Ret  3)  to  determine  reference  heating  condition:* 
and  employs  models  to  determine  the  variations  associated  with 
angle  of  attack  (tx) ,  elevon  and  body  flap  deflections  (  ,(«.  and 
JaF  )#  Reynold's  and  Mach  number  (Re),  and  sideslip  (*). 

A  program  called  HEATKST  (for  HEAT  ESTimution)  was  used 
to  reduce  the  flight  data.  The  variations  above  were  a  .mimed 
to  be  linear  derivatives  of  the  heat  rate.  Using  a  heating 
model  the  simulator  was  used  to  calculate  a  nond i me n s  1  on  a  ' 
heating  ratio.  A  thermal  model  was  also  used  within  tin- 
program  to  calculate  the  temperature  through  the  TPS.  The 
one-dimensional  differential  equations  in  this  thermal  model 
are  solved  numerically  to  propogate  the  temperature, 
sensitivity  of  the  temperature  to  each  parameter,  and  i he 
covariance  ot  the  temperature  to  the  next  discrete  time  at 
discrete  node  points  through  the  TPS  (Ret  1). 


The  flight  test  maneuvers  employ  two  techniques  to 
obtain  the  envelope  expansion  data.  A  block  diagram  of  this 
approach  is  shown  in  Fig.  2  (Ret  1).  The  first,  called 
a  Push-Over,  Pull-Up  (POPU),  is  employed  at  a  prescribed 
velocity  during  the  Orbiter  entry  when,  from  a  given  angle  of 
attack,  it  is  pitched  downward  manually  by  the  crew  at  a 
preplanned  rate,  held  at  the  minimum  angle  ot  attack  tor  five 
seconds,  then  pitched  upward  to  a  maximum  angle  ot  attack  and 
again  held  for  the  same  duration  before  being  brought  back  to 
its  original  position.  This  maneuver  is  employed  to  obtain 
aerodynamic  lift  to  drag  ratio  and  vehicle  trim  information  as 
a  function  ot  the  angle  ot  attack.  With  deviations  in  the 
flight  path  ot  opposite  signs  and  by  holding  the  maneuver  to 
only  a  few  seconds  the  original  trajectory  is  im»i 
significantly  atfected  while  at  the  same  time  AFPTC  simulator 
studies  have  indicated  that  most  thermocouple  locations  would 
be  atfected.  With  other  variables  constant,  the  heating  rate 
versus  angle  ot  attack  derivative  (Q<*)  could  be  estimated  by 
thermocouple  measurements. 

A  second  type  of  maneuver  called  a  body  flap  sweep 
consists  ot  cycling  the  body  flap  up  and  down  while  other 
variables  are  constant.  A  roll  doublet  can  also  be  performed 
while  the  flap  is  down  and  the  elevon  is  up  to  0,1.  in.  ste 
aileron  control  derivatives. 

Several  considerations  must  be  made  in  determinin',  in, 

mission  entry  profile.  A  temper  entry  could  be  p<-it  >i  

causing  higher  reference  heat  mg  but  over  a  ;.l.  <i  ter  p>  i  . 


design 


OPERATIONAL  PLIGHTS  WITH 
PLIGHT  ENHANCED 
DATA  BASE  AND  SIMULATOR 
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time  than  the  shallow  entry  that  produces  lower  peak  heat  in., 
but  higher  total  heating.  Another  coneiderat ion  is  the  r 1  up 
and  elevon  schedule.  If  the  flap  starts  at  a  low  (neat  zero) 
deflection  angle  it  should  increase  in  positive  deflection 
(downward)  to  maintain  trim.  At  the  hinge  line,  the 
previously  heated  section  of  the  flap  now  begins  to  become 
part  of  the  hinge  area  and  radiates  heat  to  the  aluminium 
chain  seal,  rib  plate,  and  other  internal  areas  in  the  a t t 
fuselage  possibly  creating  a  heating  problem.  so  the  pi.-v.-Le red 
deflection  schedule  for  the  body  flap  is  to  have  a  high 
positive  deflection  and  decrease  during  entry.  However,  a 
trade  off  must  be  made  here  with  the  elevons  which  deflect 
to  balance  the  trim.  They  also  heat  at  higher  deflections  and 
must  be  used  tor  banking  and  angle  of  attack  control. 

Execution  of  the  flight  test  maneuvers  was  successfully 
carried  out  on  the  second,  third,  and  fourth  shuttle  missions; 
STS-2,  3,  and  4.  STS-2  performed  a  series  of  maneuvers;  a 
POP U  at  Mach  20  and  three  flap  maneuvers  at  Mach  21,  17,  and 
14  as  well  as  other  aerodynamic  maneuvers.  The  hkatkst 
program  was  used  to  determine  the  heating  derivatives  at  most- 
key  locations  on  the  Orbiter.  No  problems  were  encountered  m 
estimating  derivatives  during  the  POPU  at  points  on  the  lowei 
surface  from  the  nose  to  X/L  of  .7.  The  upper  surface- 
derivatives  also  were  determined  tor  'l-  orbital  Maneuvering 
System  (OMS)  pods.  Problems  at.  a  point  on  the  elevens  and  of 
higher  than  anticipated  heating  at  one  point  on  tin-  oh;,  pet:, 
were  encountered  but  were  resolved  by  HKATKs'l  .ift.-r  som« 


manipulation. 
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Problem  and  Scope 

The  body  tlap  data,  the  heating  effects  on  it,  and  its 
attect  on  the  limits  of  the  att  Center  ot  Gravity  (CG)  <>t  tin- 
Orbiter  had  not  been  fully  reduced  until  this  report.  The 
approach  in  reducing  this  data  was  to  use  HEATKST  as 
previously  accomplished  for  other  Orbiter  locations.  This 
was  done  by  estimation  of  the  heating  ana  thermal  parameters 
tor  various  Mach  number  conditions  during  the  timeframe  when 
the  maneuvers  were  taking  place.  The  program  was  run  in 
segments  to  identify  heat  estimates  and  determine  their  trends 
throughout  the  entry.  After  these  estimates  and  their  error 
bounds  were  determined,  a  similar  procedure  was  to  be  used  on 
another  body  tlap  thermocouple  to  determine  its  effects. 
Then,  by  combining  these  results  with  results  obtained  on  the 
elevon  heating  at  the  same  time,  an  estimate  ot  the  heating 
limits  ot  the  elevon  and  body  flap  deflections  could  be 
obtained.  From  these  estimates  a  nominal  and  maximum  att  CG 
limit  tor  various  angle  ot  attack  conditions  during  entry 
could  be  produced.  This  could  establish  a  baseline  tor  limits 
(as  now  perceived)  to  the  maximum  crossrange  ot  the  orbiter 
based  on  heating  and  thermal  parameters. 

This  report  describes  the  space  Shun  le  orbiter  body 
flap  and  its  TPS  operations  during  entry.  it  also  ox,. lain 
the  theory  relevant  to  the  body  tlap  a  a  1  on  I  a  t  ion;,  i  :i  t  !.<■ 
H  IS  AT  K  ST  program.  A  detailed  analysis  ot  tin-  (<  rt  >-  ,,  o  <  : 
the  flight  test  maneuvers  conducted  th  >s  (ar  ana  the  insult  ant 
heating  on  the  body  flap  thermocouple  under  sluuy  is  in-  luu.  .. 
Results  obtained  from  those  maneuver.,  in  iol.it  a  i 
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deviation  or  the  basic  pitching  moment  trom  pre-t light  STS-1 
estimates  is  also  shown. 

Finally,  the  results  ot  the  analysis  ot  the  heating  and 
thermal  parameters  is  discussed  as  well  as  the  problems 
encountered  at  this  point  in  the  analysis.  Alternative 


approaches  and  tinal  results,  conclusions,  ana  recommendations 
complete  this  report. 


II.  The  Space.  Shuttle  body  Flap  and  tps  uy£ ration 


General  Description 

The  body  tlap  is  an  intreyal  part  ot  the  Orbiter  as  seen 
on  Fig.  3.  It  consists  of  an  aluminum  structure,  ribs, 
spars,  honeycomb,  skin  panels,  and  a  trailing  edge  assembly. 
The  main  upper  and  lower  forward  honeycomb  skin  panels  ar« 
joined  to  the  ribs,  spars,  and  honeycomb  trailing  edge  with 
structural  fasteners.  The  removable  r  forward  honeycomb 

skin  panels  complete  the  body  tlap  structure.  The  body  flap 
is  covered  with  High  Temperature  Reuseable  Surface  Insulation 
(HKSI)  and  an  articulating  pressure  and  thermal  seal  to  its 
forward  cover  area  on  the  lower  surface  of  the  body  flap  to 
block  heat  and  air  from  penetrating  to  the  structure.  The 
body  tlap  deflects  from  the  hinge  line.  physically,  the  bony 
tlap  can  deflect  +22.5  degrees  (downward)  and  -11.7  deyi 
(upward)  from  the  hinge  line.  However,  the  base  ot  the  aft 
fuselage  angles  upward  approximately  5.5  degrees  in  re  let  inn 
to  the  rest  ot  the  base  so  the  angle  the  flap  delects  is  m>i 
the  same  as  its  deflection  angle  in  relation  to  the  flow 
passing  under  the  rest  ot  the  (arbiter's  base. 

body  Flap  operation 

The  two  primary  functions  of  the  body  flap  ..r<-  '<> 

protect  the  Space  shuttle  Main  Engine's  (ss*u:'s)  it. - 1 1 :  r  y 

heating  and  to  act  in  conjunction  with  the  *•  . .  •/■.ns  as  i  •  <•  m  . •  i 
ot  gravity  trim  device.  The  arbiter  is  a  unigue  vein  i.  . « . 
that  it  is  the  first  manned  vehicle  to  reeniei  i  he  atm  ...  . 


I  FujaUpt'BotJy  Flap  Imertac* 


Fiuure.  'i  booy  t  iaj 


like  a  glider.  The  angle  ot  attack  during  entry  is  dictated 
by  the  mission  requirements.  The  Orbiter  was  designed  to  have 
a  980  nautical  mile  cross-range  maneuvering  capability.  Since 
DOD  planned  launches  trom  the  Western  Test  Range  at  Vandenberg 
ARB  presently  need  this  capability  to  fulfill  their  nominal 
and  Abort  Once  Around  requirements,  the  Orbiter  must  enter  at 
an  angle  of  attack  of  36  to  38  degrees  and  roll  angle  ot  60 
degrees.  Later  in  the  entry,  to  obtain  more  cross-range  by 
increasing  the  lift  to  drag  ratio,  the  angle  of  attack  is 
reduced  to  as  low  as  26  degrees.  Unfortunately,  this  also 
severely  increases  the  heating  at  many  locations  including  the 
tlap.  This  is  due  to  higher  Reynold's  number  and  body  flap 
deflection  at  lower  angle  ot  attack.  The  Map  deflects  more  to 
maintain  trim  at  lower  angles  ot  attack.  The  Reynold's  number 
increase  is  due  to  a  higher  velocity  and  density  with  the  lift 
coeticient  rising  taster  than  the  drag  coeticient  at  a  lower 
altitude  in  the  flight  corridor.  Later  in  the  entry  a  gradual 
pitch  down  is  initiated  starting  below  Mach  m. 

The  elevons  and  body  flap  also  must  be  used  together  to 
control  the  trim  of  the  Orbiter  during  the  entry  period.  it* 
nominal  CG  position  is  .667  trom  the  nose  tip  while  the  most 
aft  CG  position  is  .675.  The  aft  CG  increases  body  t  I 
deflection  ( I  )  and  elevon  det  lection  (  )  tor  trim.  it  is 

therefore  important  to  know  the  heating  rates  on  the  body  f  l  .ip 
and  elevons  at  various  deflections  to  determine  i:  t  t . » 
reuseable  surface  insulation  on  both  can  protect  them  and  I  < 
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reused 


The  Thermal  Protection  System  (Tl-s)  on  the  Orbit  >.-r  w a  ■; 
designed  to  protect  the  Orbiter  and  still  be  reuseable  rm  .a 
least  one  hundred  missions.  The  TPS  on  the  Orbiter  consists 
primarily  ot  four  types  of  materials,  each  designed  to 
insulate  the  Orbiter's  aluminum  skin.  The  High  Temperature 
Reuseable  Surface  Insulation  (HRSI)  is  able  to  withstand 
temperatures  as  high  as  2b00  R  tor  13 ne  mission  and  be 
reuseable  when  not  exposed  to  temperatures  higher  than  2270  R. 
The  tile  consists  of  amorphous  borosilicate  glass.  To  prevent 
damage  to  the  fragile  tile,  strips  ot  felt  padding  ate  used 
between  the  tiles  and  the  orbiter's  skin.  This  strain 
Isolation  Pad  (SIP)  is  bonded  to  both  the  tile  and  the  orbiter 
with  a  red,  room-temperature  vulcanizing  ( RTV )  cement  with  an 
iron  oxide  base.  Other  strips  of  padding  called  fillet  bars 
are  placed  between  each  tile. 


A  program  named  HEATEST  (tor  HEAT  ESTimating)  was  used 
to  reduce  flight  test  data.  A  flow  chart  of  HEATEST  is 
depicted  in  Fig.  4.  HEATEST  incorporates  systems  theory  of 
theoretical  models  in  order  to  compute  parameters  that  effect 
the  heating.  The  theory  from  scientific  method  has 
established  a  linearly  independent  set  of  parameters.  Due  to 
the  number  of  these  parameters,  the  systems  theory  in  HEATEST 
combines  some  of  these  variations  in  order  to  compute  them 
quickly.  The  variations  are  assumed  to  be  linear  derivatives 
of  the  heat  rate.  The  heating  ratio  Q/Qref  was  assumed  to  be 

Q/Qref  =.  f  («*,  *  ,<}  ,**)  (1) 

where  Qref  is  the  dimensionless  reference  heat  rate  on  a  one 
foot  diameter  sphere  given  by 

Qref  *  17700^  V^/loV”  (1-Hw/Ho)  (2) 

hw  *  .24  [Oref/(<re)J  M  (3) 

ho  *  .24  T„+V,//5006  3  (4) 

where  the  heating  rate  is  in  British  Thermal  units  (BTU's)  per 
toot  squared  per  second,  r  is  the  Stef an-Bol tsman  constant 
(4.761x10''*  ),  £  is  the  emissivity,  /J.  is  the  atmospheric 
density  in  slugs  per  cubic  foot,  is  the  relative  velocity 
in  teet  »r  second,  and  is  the  atmospheric  temperature  in 
degrees  Rankine. 
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Kicjure  4.  ied  HfcATKST  Klow  Diagram 


This  heating  model  and  wind  tunnel  data  at  specific 
Orbiter  TPS  points  were  used  to  calculate  the  simulator 
heating  rate  ratio  using  linear  interpolation  routines  tor 
real-time  processing  ot  aerodynamic  ciata.  For  the  HEATEST 
data  reduction  program  the  heating  ratio  took  the  form 

0  =  Q/0 ref  -  0o  +  0^ (« -  *0)  +  0,(P  -».) 

+  Qiv  -  i-F,  )  +  (U-  U.) 

+  0  |  [  log  (  Re )  -  log  (  Re„  )  ] 

♦  QM(M*-M»a)  (5) 

where  Qo  is  the  magnitude  at  the  reference  conditions 
specified  by  the  variables  with  the  zero  subscripts.  In  the 
body  flap  calculations  the  &  terms  and  M  terms  were  neglected. 
The  subscripts  on  the  heating  ratio  (0)  represent  partial 
derivatives  with  respect  to  the  variables  ot  angle  ot  attack 
(a<),  sideslip  (0),  the  logarithm  to  the  base  ten  ot  the 
Reynold's  number  (Re)  based  on  its  characteristic  length, 
elevon  deflection  angle  (Sc)i  body  flap  detlection  angle 
and  Mach  number  (M«.). 

These  heating  derivatives  are  essential  in  the 
evaluation  ot  determining  what  factors  influence  the  entry 
heating  at  the  flap.  Another  value  ot  Qret  could  be 

-  S  a/3 

Qret  -  .  33 2(  To-Tw )  Re  '  V^V-Cp/i  778 . 3Pr  )  ( «> ) 

which  is  based  on  Eckert  theory  for  laminar  clow  with  no 
deflection.  Laminar  flow  is  approprate  tor  most  areas  on  the 
Orbiter.  However,  as  will  be  later  demonstrated,  it  is  not 
correct  tor  the  flap  and  other  regions  where  the  turuui  >nt 


equation  for  the  deflection 


Qref  »  .0296(To-Tw)Re“  ^0V-Cp(T*/To)/(778.3Pra/,)(P/P(O)  (7) 

and  the  T*  value  is 

T*  =  Ta  +  .5(Tw-T,,  )+.22(TR-Ta  )  (8) 

is  more  representative  of  the  heating.  TR  is  assumed  as  .9To. 

The  change  in  the  heat  rate  (AQ)  from  the  reference 
conditions  is  given  by 

AQ  =  0  -  Qo  (9) 

where  the  magnitude  of  the  reference  condition  is  subtracted. 
The  assumption  is  made  that  for  short  time  segments  these 
derivatives  are  constant  and  the  lateral  conduction  in  each 
tile  is  small. 

Thermal  Mode-L 

The  AFFTC  simulator  also  employs  a  thermal  model  ot  the 
TPS  tile  to  calculate  the  temperature  at  various  locations 
within  the  tile  and  its  surroundings.  A  description  ot  this 
model  is  shown  in  Fig.  b.  The  interior  ot  the  tile  is  divided 
into  equal  thickness  and  the  entire  structure  ot  different 
materials  is  designated  by  blocks.  Numerical  solution  ot  the 
resulting  equations  resulted  in  an  accurate  simulation  ot 
surface  and  bondline  temperatures  at  specific  locations.  ihi- 
honeycomb  also  is  included  in  this  analysis.  There  are  two 
face  sheets  on  the  honeyccnb  with  vecy  little  aluminum  in 
between  so  that  there  is  radiative  heat  transtet  caused  oy  tin. 
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Kiyure  5.  THb  Model  Cross  .section 


different  temperatures  between  the  two  sheets.  it  it  is 
assumed  that  there  is  no  material  present,  then  in  it  is  only 
the  radiation  between  the  two  node  points  or 


(l-AKcrcr,*  -t./) 


( 10) 


where  A  is  area  per  unit  area  aluminum  between  the  tacesheets. 

In  the  tile,  node  points  are  spaced  equally  into  sma  L  L 
elements  of  length  AX  for  a  total  Oa  t,  node  points.  Each 
block  has  four  equally  spaced  nodes.  Each  block  represents  a 
different  material  with  different  properties.  Although 
thermocouples  can  be  placed  at  various  node  points  in  each 
block,  for  our  purposes  the  only  thermocouple  is  located  at 
the  second  node  (i*2)  which  is  AXa  from  the  surface.  The 
convective  heat  rate  (0)  input  at  node  1  has  most  of  the  heat 
radiated  away  and  conducts  a  small  amount  to  node  2  and  the 
rest  of  the  tile.  Block  C  is  bonded  to  the  tile  by  kl'V 
adhesive  to  the  SIP  which  is  also  attached  to  the  honeycomb 
face  sheet  by  RTV. 

An  ordinary  differential  equation  tor  the  temperature 
(Ti)  at  the  ith  node  point  was  obtained  from  an  energy  balance 
tor  each  element.  A  system  of  L  nonlinear  differential 
equations  results  (Ret  1). 

The  energy  balances  tor  each  element  are  summed  together 
to  form  the  heat  transfer  equations  (Ref  4) 


1-1  <C,V>V2»T  *  -K«''4XAT,+l‘./4XA'h 

-rCIT,''  -T-h  +  f  (ec.P.4  ,#e)Qret 


i*2  (Ca/aAXa/2  +  Cj^AXa/2)Ta-  Ka/AXaT3  (12) 

-(Ka/AXa  +  KA/AXA)Ta  +  K#/AX,T3 

i  Ca/#  (AX  +AX  j)/2T{  -  K./AX,.,  T,_,  (LI) 

-(Kf/AXj./  +  Kb/AX  j  )T  (  +  Kg/AXj  Tj., 

i-L-2  (C,  »%AXl.3  /2  +  Ctf/cAXc/2)f  ,.a  =  Kc/  AXCTL.3  (14) 

- ( Ka/AX  L.j  +  Kc/AXc)Tl.2  +  KC/AXCT 

i”L-l  ((-sir  ^si*AXtlp/2  +  CjfjV-fnr vAX^rv  (IS) 

+C  HC -/KCAXHCA/2  )i  , 

-K0/(AXc)TL.Jl  “( Kp*  +K0‘)Tt.(  +KT  L 

+  CMco-<  [Tt+460]  -  (T  ,  +46U)  ^  )  ( 1  -  A  ) 

i=L  (CMC/HCAXrACt+CH£^cA  XMC/2  (lb) 

+cDy„(Axa/2  +AXa  ,  N  k)t  L 
(Tt+460)4'-(TO+460)  <,)-HF(Tt  +  460-TeO)+KD"TL.  , 
-KBTL+£Mccrl(TL+460)^-(TL.(  +460  )*  ( 1-A)  ) 

The  thermal  conductivity  ( K )  is  a  function  of  temperature  and 
local  pressure.  The  conductivity  is  evaluated  at  an  average 
temperature  between  the  nodes.  The  specific  heat  (C)  is  also 
dependent  on  the  temperature  at  the  node  point.  The  dot 
represents  the  time  rate  of  change  of  the  temperature  (Ti)  at 
each  node  with  a  distance  (AXi)  between  each  node  and  is 
approximated  by  a  first  order  backward  difference. 

If  the  surface  conditions  can  be  specified,  i.e.  the 
thermal  emissivity  (€),  the  coating  thickness  (Axa),  .mu  the 
aerodynamic  heating  rate  t(<*,  0  ,  &  ,Re) ,  the  temperature  through 
the  tile  can  be  calculated  by  solving  the  system  of  ordinu, 


differential  equations.  The  specific  heat  and  density  (/”)  are 
assumed  to  be  known  from  the  design  but  emissivity  can  vary 
from  .75  to  .92  and  the  thermal  conductivity  (a  function  of 
tile  pressure)  is  variable  with  flight  pressures  inside  the 
tile  in  doubt  (Ref  1  and  4).  Using  a  stochastic  estimation 
process  the  unknown  parameters  are  found  from  systems 
identification  theory. 

All  temperatures  at  each  node  through  the  TPS  are  not 
measured  in  the  previous  described  model  because  thermocouples 
are  not  located  at  each  node  point.  Thus  an  intermediate  step 
uses  an  extended  Kalman  filter  to  estimate  all  node 
temperatures  whenever  a  thermocouple  sample  was  available. 
This  was  based  on  the  thermocouple  measurement  and  heating, 
thermal,  and  error  models.  This  filter  also  filtered  out 
noise.  The  parameters  were  then  updated  by  a  gradient 
algorithim  to  maximize  a  maximum  likelihood  function  tor  each 
parameter.  one  or  several  of  these  parameters  were  estimated 
over  a  planned  time  segment.  Thus,  it  the  segment  is  very 
long,  a  nonlinearity  in  these  parameters  over  the  entire  entiy 
can  occur.  Parameters  are  used  to  update  the  simulator  and 
therefore  data  is  enhanced  by  flight  data. 


One  method  of  approximation  of  the  initial  temperature 
Ti  is  based  on  an  empirically  determined  RC  time  constant  an*', 
on  the  radiation  equilibrium  assumption.  The  circuit  analogy 
is  used  to  compute  an  equilibrium  temperature  (Teq)  of 


Teq(tn)  ■  gT, < tn)-(g-l)T, ( tn- •) 


(IV) 


where 


g  -  1/ {l-exp(-At/RC) ) 


( IB) 


and  the  RC  time  constant  value  can  be  equated  to  the  lump 
parameter 

RC  -  CpAX  */K  (  19  ) 

This  varies  with  the  location  in  the  TPS.  In  the  coating  kc 
is  about  1/4  second  while  empirically  RC  for  the  tile  is  about 
3  seconds.  The  equilibrium  temperature  works  well  in  t.n«* 
program  when  you  take  0  and  assume  equilibrium.  The  technique 
of  using  the  RC  analogy  is  very  crude  and,  when  first  derived, 
was  not  intended  for  use  in  flight  data  correlation.  In  tins 
report  it  is  being  used  in  HEATEST  as  an  approximate  method  ot 
obtaining  initial  conditions,  establishing  trends,  and  quick 
turnarounds. 

In  this  approximation  program,  conduction  effects  are 
approximated  by  the  RC  circuit  analogy.  It  the  temperature 
changes,  the  RC  approximation  amplifies  the  heating  value  but 
is  within  reasonable  accuracy.  The  RC  approximation  ampLities 
spikes  (noise  caused  by  the  8  bit  words)  and  only 
approximates  transients  such  as  angle  ot  attack  effects  when 
using  RC* 1 .  The  RC  network  analogy  is  explained  in  Appendix 
2 . 


21 


To  expand  the  Orbiter's  envelope  and  remove  placards, 
flight  test  maneuvers  were  planned  for  STS-2  and  4,  The  AFFTC 
used  a  systems  approach  to  the  envelope  expansion.  The 
placards  related  to  the  body  flap  included  heating  of  the  flap 
and  elevons  at  various  Mach  numbers  and  angles  of  attack, 
elevator  and  body  flap  effectiveness,  verification  of  the 
basic  pitching  moment  curve  (Cmo),  and  lift  to  drag  ratio  in 
the  hypersonic  regime. 

To  do  this  a  POPU  maneuver  from  45  to  35  degrees  angle 
of  attack  was  planned  on  STS-2  at  Mach  20  and  three  body  t Lap 
sweeps  at  Mach  21,  18,  and  14.  STS-4  planned  a  PUPO  maneuver 
at  Mach  14  and  at  Mach  8.  The  flap  was  cycled  through  sweeps 
from  17.1  to  7  degrees  positive  (downward)  deflection. 

Performance  of  Maneuvers  on  STS-2  and  4. 

The  entry  profiles  of  STS-2  and  4  occureJ  as  planned. 
Thermocouple  data  was  recorded  on  board  and  telemetry  sent  in 
real-time  when  in  view  of  a  station.  Unfortunately,  critical 
data  was  lost  on  STS-4  when  the  recorder  tailed  and  real-time 
data  was  not  acquired  until  after  most  of  the  high  Mach 
heating.  Data  from  the  body  flap  thermocouple  used  in  t  lie- 
data  reduction  tor  this  report  is  shown  in  Fig.  6.  The  STS- t 
and  4  thermocouple  data  at  the  same  location  is  also  shown  !<»r 
comparison.  The  elevon  schedule  tor  the  first  tour  missions 
has  been  ie*-l,  1,  3,  and  5  oegrues  respect  i ve  1  y  with  u  umn. 
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expanding  high  Mach  number.  At  the  same  point  the  6  er  has 
been  16,  14,  12,  and  4  degrees  respectively,  only  the  STS-2 
data  is  complete  enough  to  do  a  thorough  analysis.  The  t  lap 
and  elevon  time  histories  during  the  STS-2  maneuvers  arc- 
displayed  in  Fig.  7.  A  table  of  STS-2  maneuver  events  is 
shown  in  Table  1.  The  heating  at  the  body  flap  thermocouple 
location  shown  indicated  a  significant  increase  in  temperature 
over  that  anticipated.  It  should  be  noted  that  this  was  not 
the  only  location  where  anticipated  thermocouple  temperatures 
varied  from  actual  results.  Some  locations  displayed 
significantly  lower  temperatures. 

The  orbiter  flight  control  software  logic  contains  a 
body  flap-elevator  interconnect  which  is  intended  to  maintain 
the  elevator  on  a  preplanned  schedule  as  a  tunction  ot  Mach 
number.  On  all  tour  test  flights  it  automatically  trimmed  the 
body  flap  at  an  angle  much  higher  than  anticipated.  it  was 
concluded  that  the  error  was  in  the  basic  pitching  curve  t’mo. 
This  information  was  reduced  (Ret  5)  to  a  series  of  curves  us 
shown  in  Fig.  8  and  9  comparing  the  maximum  elevon  positive 
deflection  (in  order  to  retain  its  other  requirement  of  ent i / 
pitch  control)  and  body  flap  position  for  various  Mach  numbers 
and  angles  ot  attack  at  nominal  and  maximum  ate  CG  positions. 
This  data  indicates  that  it  is  important  to  know  the  therm.ii 
limitations  on  the  flap  and  elevons  to  determine  it  mission 
requirements  can  be  met.  other  data  obtained  conclude.;  Umi 
the  body  tlap  and  elevon  effectiveness  as  well  as  hypei  son  i  ■  - 
lift-to  drag  ratio  were  the  same  as  predicted  (Ret  b). 


Table  1  STS-2  Flight  Test  Maneuver  Key  Events 


EVENT  MISSION 

TIME 
(SEC  OF 
OMT) 

TIME 
FROM 
k  £ 

(SEC) 

ANGLE 

OF 

ATTACK 
(  DEG) 

LOG  Re 

THERMOGi 
TEMPERA 
(  DEGHKI 

Entry 

Interface 

75037 

00. 

41.05 

2.365 

117 

1st  Flap 
Sweep  (+) 

75742 

705. 

39.86 

6.089 

1844 

Negative ' 
Deflection 

75786 

749. 

38.32 

6.146 

2054 

End  Sweep 

75812 

775. 

38.12 

n  .  14  6 

1713 

POPU  (PO) 

75830 

79  3. 

38.91 

6.166 

20.13 

Minimum 

75845 

808. 

34.59 

6.178 

1992 

POPU  (PU) 

75846 

809. 

34.59 

6.179 

1992 

Maximum 

75865 

828. 

45.59 

6.200 

2197 

End  POPU 

75887 

850. 

38.58 

6.220 

2U65 

2nd  Flap 
Sweep  (-*•) 

75994 

957. 

40.21 

6.375 

2197 

Max.  + 
Deflection 

75996 

959. 

39.59 

6.38 

2228 

Negative 

Deflection 

75997 

960. 

39.60 

6.382 

2258 

Peak 

Heating 

75998 

961  . 

39.40 

6.383 

2278 

Maximum  76008 

-  Deflection 

970. 

4  1.41 

6.395 

223H 

End  Sweep 

76015 

978  . 

41  .00 

6.401 

2126 

3rd  Flap 
Sweep 

76128 

lO^l  . 

19.81 

6.6  38 

21  3«> 

Max .  4- 
Deflection 

76138 

110  1. 

38  .93 

b  .  6 1.  ' 

215' 

Max .  - 

Det lect ion 

76172 

113  5. 

39 .60 

b  .  7  3  1 

20?  1 

End  Sweep  76181 
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The  fixed  body  flap  during  the  POPU  made  it  possible  to 
identity  the  0^  and  Qo  using  HEATEST  it  the  0  *«,  derivative  is 
known.  Other  selected  parameters  can  be  expressed  to  be 

QO,Q<x/Q  'Oi  l  Xs  (^(J) 

LOG  ke  '  i,F 

where  AXa  is  the  effective  thermocouple  depth  or  coating 
thickness.  The  goal  of  the  data  reduction  program  in  Hfc.AThSi 
is  to  obtain  best  estimates  of  these  parameters  during  the 


maneuver 


V.  Results 


Initial  Computer  Runs  and  Analysis 

The  STS-2  flight  thermocouple  data  tor  the  Location  on 
the  body  flap  shown  earlier  was  input  into  the  HEATEST 
program.  Also  input  was  a  best-Est imated-Tra jectory  (BET)  and 
atmosphere  from  NASA  Langley  kesearch  Center  tor  reducing 
"model  error"  caused  by  bias  in  the  heating  ratio  from 
atmospheric  density  error.  A  problem  of  time  skews  in  th< 
data  for  other  points  was  not  apparent  in  the  flight  data 
input  at  this  location.  There  was  however  some  drop  outs  of 
the  body  flap  deflection  angle  data  during  one  body  flap  sweep 
that  had  been  "tilled  in"  prior  to  its  use  with  this  study, 
both  of  the  earlier  described  errors  can  produce  errors  in 
HEATEST  estimates  which  cannot  be  identified  by  HEATEST  as 
such.  The  procedure  developed  for  running  the  proyrain  was  to 
take  the  entry  data  and  to  break  up  the  computer  runs  into 
segments.  Each  segment  covered  a  particular  region  of  the 
entry  where  a  maneuver  was  occuring.  These  segments  w«-n- 
further  broken  down  into  smaller  time  segments.  For  example, 
a  run  called  HEAT  3  was  computed  between  STS-2  times  ot  757  u 
and  75811  seconds  GMT,  a  segment  after  entry  interface  (Ell  <ii 
400,000  ft.  altitude.  This  is  where  the  first  body  tl  1 1  S  W  « *  ♦  •  i  1 
took  place.  This  was  further  segmented  in  attempt,  to  h.ivi- 
HEATEST  estimate  individual  parameters  and  uncertainty  :  nine  . 
These  best  estimate  derivatives  or  slopes  were  then  used  in 
later  time  segments  to  determine  it  the  model  would  "match" 
the  variations  in  temperature. 


From  experience  gained  by  others  (Ret  6),  derivatives 

which  had  the  greatest  impact  were  0  f,  Qo,  Q ■  ,  Q i  ,  and 

LOG  Re 

the  smallest  Qj,  .  The  coating  thickness  (AXa)  was  initially 
assumed  to  be  .002  feet  and  later  confirmed.  For  most  runs, 
the  temperature  was  calculated  to  obtain  an  approximate 
initial  condition  of  the  surface  te  -perature  usiny  an  RC 
circuit  analogy.  A  best  guess  of  all  derivatives  at  an  initial 
time  segment  was  also  made  to  improve  initial  conditions.  The 
next  segment  of  time  was  taken  ideally  when  only  one  of  the 
heating  parameters  (ex,  l* )  was  changing.  At  this  point 

the  program  was  directed  to  allow  the  derivative  of  the 
changing  parameter  to  "float"  to  a  new  value  estimated  by  the 
program.  This  value  was  then  "locked  up"  or  fixed  during  the 
next  time  segment  in  order  to  establish  the  parameter  values 
in  that  region.  As  this  was  taking  p ’ ace  the  Kalman  filter 
was  establishing  a  new  surface  temperature  based  on  its 
updated  model  parameters.  It  the  parameters  were  correct,  t  ti«; 
apriori  and  aposterori  temperatures  would  match  the  variations 
in  the  flight  thermocouple  temperature  each  second  within  an 
uncertainty  of  approximately  ten  degrees,  which  is  the 
accuracy  of  the  8-bit  word  measurement  resolution.  Although 
the  &e  parameter  is  significant  at  other  Orbit  or  flap  locations 
near  the  centerline,  its  value  here  is  very  small  and  was 
given  a  small  fixed  value  during  most  runs.  It  was  also  done 
in  an  attempt  to  prevent  trading  of  its  value  with  othei 
parameters.  Otherwise,  with  the  continued  variation  o*  tin 
elevon  during  entry,  the  <j^e  would  have  been  required  to 
"tloating"  through  most  of  the  maneuver  times  studied. 


Using  the  previous  technique  was  for  most  cases 
unsuccessful.  The  temperatures  would  drift  as  more  segments 
were  taken  within  each  run  and  a  good  match  tor  a  lomj  time 
segment  was  never  truly  established.  Trading  of  derivative 
values  such  as  Qo  and  or  QjBf  and  would  occur  it  more 
than  one  derivative  was  being  estimated  at  the  same  time. 
Since  both  parameters  may  be  changing  at  the  same  time  then 
both  would  require  estimating.  It  the  values  derived  from  the 
program  were  the  correct  ones  the  trading  may  still  occur. 
Examples  of  the  parameter  values  obtained  and  t  fie  re 
uncertainty  bound  during  different  portions  ot  the  entry  are 
shown  in  Fig.  10.  The  first  graph,  depicting  a  search  tor  the 
right  Qp<  derivative,  were  all  taken  at  various  portions  ot  the 
POPU.  The  second  graph  depicts  a  search  tor  the  proper 
derivative  and  were  all  (except  the  13  degree  run)  done  during 
the  three  body  flap  sweeps.  The  third  chart  again  looks  tor 
the  right  Q«  during  various  portions  in  the  entry.  Normally 
values  with  large  error  bounds  are  not  used. 

The  results  of  the  attempts  to  calculate  angle  ot  attack 
and  body  flap  deflection  effects  are  not  conclusive  but  those 
with  a  small  error  bound  appear  to  have  approximately  the  same 
values  which  may  indicate  that  these  are  the  correct  parameter 
values.  Only  the  trend  in  the  change  ot  these  values  (it  'hey 
do  change)  during  the  entry  may  be  all  that  is  required  to  oo 
once  HEATRST  can  be  successfully  run  through  the  maneuvers. 

It  was  assumed  at  this  point  tha*  t'.*  problem  was  caused 
by  one  or  two  items;  the  Reynold's  number  derivative  and/or  a 
nonlinearity  in  the  other  derivatives,  or  flap  flow  effects. 


The  Reynold's  number  derivative  is  unique  in  that  t  low 
transition  causes  a  rapid  increase  in  its  generally  small 
value  and  then  decreases  rapidly  when  tully  turbulent.  This 
effect  can  also  cause  difficulty  in  interpreting  the  results 
if  the  transition  occurs  during  the  maneuver.  If  interpreted 
correctly  it  indicates  sensitivity  to  transition  in  one  ot  the 
variables,  particularly  it  the  flow  returns  to  a  laminar 
state.  Such  a  condition  occured  on  STS-4  on  the  lower  surrace 
centerline.  The  assumed  Reynold's  number  effect  could  also  be 
a  Mach  number  effect  and  real  gas  effect  but  the  two  cannot  be 
separated  during  a  given  flight  other  than  by  theory.  The 
Mach  number  effect  could  also  be  pulled  out  in  lieu  of  ke. 
The  Reynold's  number  changes  so  slowly  that  the  only  way  to 
accurately  determine  its  value  is  to  run  the  program  for  as 
long  as  possible  (over  100  seconds  is  best)  with  no  other 
changes  taking  place  in  the  region.  U  fortuately,  the  region 
in  question  is  very  dynamic  and  the  program  as  established 
could  not  distinguish  Reynold's  number,  Qo,  and  effects  so 
the  slope  could  not  be  obtained.  When  the  Reynold's  number 
runs  were  made  these  too  would  not  converge  except  at  one 
point  in  the  entry.  This  single  point  was  not  sufficient  to 
establish  a  trend,  figure  11  depicts  tike  0n>,  vetse  i.oti  ko 
results  obtained.  Although  the  error  bounds  appear  small  it 
must  be  remembered  that  derivatives  are  trading  value.-, 
resulting  in  a  small  error  band  tor  an  incorrect  value. 
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l-q£  Ha  a  tiny  Models 

w  Two  methods  of  removing  nonlinearity  effects  are  to 

change  time  segments  in  HEATEST  or  change  the  model.  Many 
time  segments  were  tested  without  success.  In  systems 
identification  theory  model  error  is  a  major  concern  and 
problem  (Ref  1).  Since  the  model  had  worked  with  very  little 
difficulty  in  pther  Orbiter  regions  it  was  surprising  and 
disturbing  that  it  occured  at  this  point.  The  only  option  at 
this  point  was  to  investigate  other  more  complex  models  based 
on  more  complex  theories  and  compare  them.  In  comparing  these 
theories  it  was  not  the  intent  to  say  they  were  totally 
accurate  and  therefore  could  be  used  in  predicting.  The  main 
objective  was  to  use  various  theories  to  see  it  they  agree 
with  the  flight  data  in  magnitude,  slope,  or  linearity. 

1 

In  first  plotting  the  flight  ccnd.'ion  Q/<jce£  versus  Log 
Re,  obtained  by  only  estimating  Qo  in  HEATEST  with  assumed 
I  value  of  0«k  and  0^,,  ,  the  Reynold's  analogy  theory  was 

compared  for  laminar  and  turbulent  tlow.  This  is  shown  in 
Eiy.  12.  Hand  calculations  revealed  that  the  turbulent  tlow 
I  more  closely  matched  the  flight,  data.  Since  the  theory  in  the 

turbulent  case  uses  Re’*  an  axis  change  of  o/Orei  ver  se  ke  '* 
revealed  a  more  linear  slope.  Calculations  of  these  result:, 
and  others  are  in  Appendix  A. 

Calculations  using  the  normal  shock  tables  and  keynote's 
Analogy  as  Rockwell  established  in  their  data  book  (her  /) 
were  also  compared.  These  values  did  not  match  in  magnitude 
or  slope.  In  the  region  of  interest,  the  magnitude  w.,s  t  > ». . 


high  at  the  start  and  too  low  at  the  end. 


>/ Orel-  For  Flight  Data  ana  Reynold's  Analogy 


In  another  approach,  it  is  considered  we  empirically 
know  that  down  the  centerline  of  the  Orbiter  prior  to  t  li<-; 
boattail  expansion  that  the  flow  is  approximatel  y  Mach  i  (i<er 
6).  The  stagnation  and  after  the  shock  tempe r a tu res  can  he 
evaluated.  From  these  values,  the  conditions  at  the  flap  and 
more  accurate  values  of  viscosity,  Reynold's  number,  ami 

O/Qret  can  be  calculated.  The  flap  angle  can  be  fixed  verses 

8  8 

changes  in  time  and  Re ‘  .  At  a  given  Re’  the  flap  can  be 
changed  in  deflection  angle  and  Q/Qret  changes  due  to  this 
deflection  can  be  determined.  The  minor  changes  in  the  angle 
of  attack  are  assumed  to  be  linear  and  the  elevon  changes  are 
not  significant  at  this  point  on  the  flap.  The  calculated 
results  again  did  not  match  in  slope  but  were  close  to  the 
flight  data  in  magnitude. 

It  was  stated  in  a  recent  AIAA  paper  (Ref  8)  that 
Newtonian  flow  analysis  could  be  used  for  the  lower  surface  of 
the  Orbiter  in  this  hypersonic  regime.  Calculations  using  t h i s 
approach  concentrated  on  determining  the  specific  heat  at  high 
Mach  numbers  (Ref  9)  and  low  density.  The  flap  angle,  weak 
shock  expansion,  and  shock  tables  determined  the  downstream 
Mach  number,  density,  velocity,  Reynold's  number,  and 
viscosity.  When  calculating  from  these  values,  Q/oret  tor  We 
values  were  not  close  in  slope  to  the  flight  data.  However, 
the  magnitude  of  values  was  close  to  the  flight  results.  It 
was  observed,  in  this  analysis  and  the  pi ev ions  weak  shook 
calculations,  that  variations  in  body  t  lap  at  a  given  Ha«.h 
number  and  angle  of  attack  did  not  change  significant  !y  t  he 
O/yret  values.  The  Newtonian  analysis  however  rev.-,,le-i  tii.u 


the  changes  in  angle  of  attack  did  significantly  influence 
the  slope  of  the  O/Qref  values. 

A  fifth  approach  was  to  consider  the  Eckert  flat  plate 
formulas  (Ret  10)  used  in  the  program  to  calculate  the  heating 
ratio  Q/Qret  but  now  adding  the  downstream  conditions  and  the 
deflection  angle  terms.  Using  the  turoulent  flow  equation 

,  f 

and  plotting  the  resulting  y/Qref  values  verses  ke  tor  the 
same  entry  conditions  as  previous  calculations  resulted  in  a 
series  of  points  slightly  higher  than  the  flight  data  at  the 
beginning,  crossing  the  flight  values  during  the  first  body 
flap  sweep  and  POPU  maneuvers,  then  leveling  oft  at  values 
below  the  flight  data  later  in  the  entry. 

The  sixth  approach  was  to  consider  another  model  that 
would  use  the  approximate  downstream  conditions  in  calculating 
the  heat  rate  and  remove  the  influence  of  Reynold's  number 
completely  from  the  problem,  by  knowing  that  the  Q  value  in 
the  Newtonian  and  Reynold's  analogy  is  a  function  of  Stanton 
number  and  enthalpy,  as  well  as  the  Reynold's  number  and 
distance  from  the  stagnation  point,  the  Qret  equation  can  be 
rearranged  to  be  a  function  only  of  the  downstream  density, 
velocity,  and  enthalpy  with  the  distance  and  Stanton  number  a 
constant  (Ret  11).  Enthalpy  could  be  and  was  calculated  using 
the  freestream  velocity  or  the  specific  heat  and  stagnat  ’.«:i 
temperature  tor  the  input  conditions  used  previously.  lit¬ 
re  s  u  1 1  s  were  encouraging.  Compared  with  the  i  .  t  !  i  - i 
approaches  and  flight  date.,  the  results  followed  in  c.ws* 
agreement  with  the  flight  da*  a  through  the  first  two  hi.um*uvi  •> 
and  its  magnitude  matched  t  lie  flight  data  y/yre  r  at  t  r  ■*  i  tit  i 


body  flap  sweep.  It  roust  be  remembered  that  the  flap 
detlection  tor  these  calculations  is  considered  constant  and 
therefore  no  body  flap  sweeps  are  input  into  the  calculations. 

The  seventh  and  final  method  to  be  examined  was  m 
approach  used  by  NASA  in  their  analysis  of  the  heating,  t  hi- 
Spalding  and  Chi  technique  (Ref  12).  This  method  finds  t  In- 
value  of  the  skin  friction  coeticient,  Ct.  To  use  this 

technique,  the  flap  conditions  were  assumed  to  be  as  they  wer«- 
in  the  weak  shock  calculations.  Assuming  that  the  flap  prior 
to  the  weak  shock  experiences  a  tree  stream  Mach  number  ot  j, 
the  edge  condition  experienced  by  the  flap  becomes  Ma  =  1.27  tor 
a  flap  detlection  ot  +7  degrees.  The  procedure  described  in 
the  reference  was  performed  tor  the  seven  entry  points 
previously  used.  The  result  was  an  almost  linear  reduction  ot 
the  Ct  value  from  .99  to  ,b0  from  that  ot  the  approximation  ot 
.006  used  in  the  Stanton  number-Enthalpy  calculations.  When 
this  value  was  applied  to  those  calct’.d  ions,  the  resulting 
Q/Qret  magnitude  was  almost  identical  to  the  flight  data 
before,  during,  and  through  the  POPU  maneuver.  However, 
because  the  flap  deflections  were  not  considered,  the 
resulting  values  did  not  reflect  a  change  during  the  two  bogy 
flap  sweeps. 

Figure  13  graphically  compares  all  '>t  these  methods.  The 
Reynold’s  analogy,  Eckert,  and  Kntha lpy-Stanf on  number  values 
with  and  without  the  Spalding  and  Chi  Ct  values  at<-  m<n< 
closely  compared  with  the  t  light  data  in  Fig.  14.  Kig  ii.  is 
displays  the  wind  tunnel  data,  the  NASA  simplified  i  ■  ■ » *  1  , 

model  and  the  weak  shock  results  using  M  -  i  prior  to  t  i •  t  !  i . 
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Q/Oret 


The  Enthalpy-Stanton  number  approximation  combined  with 
the  Spalding  and  Chi  Ct  values  appears  to  be  closest  io  i  nc 
flight  data  but  would  be  extremely  ditl  icuIL  to  use  in  the 
HKATEST  program.  The  primary  difficulty  is  in  determining  in*: 
downstream  conditions  accurately  at  the  flap.  This  is  also 
true  ot  the  other  approaches  attempted  except  tor  the 
Reynold's  analogy  and  the  Hokert  equation  not  using  the 
deflection  correction  factor  and  the  temperature  ratio  t*/ to. 
It  is  obvious  that  the  weak  shock  and  Newtonian  flow  anjLyst.. 
have  a  similar  problem.  Placing  the  new  Ct  values  in  in. 
Newtonian  and  Weak  shock  calculations  made  no  sign  it  leant 
change  in  their  lack  ot  accuracy  compared  to  the  flight  data. 
Invalid  assumptions  made  or  calculations  performed  altered 
the  downstream  Reynold's  number  in  the  weak  shock  ca  leu  1  at  i  on-- 
and  the  magnitude  ot  the  Q/Ore t  in  the  Newtonian  analysis. 

In  the  calculations  tor  the  downstream  conditions  sever  a  ■ 
assumptions  must  be  made.  First,  that  the  treestream  ' low 
past  the  expansion  at  the  boattail  reaching  the  flap  is 
approximately  Mach  3.  Second,  that  the  ,  :■  w  is  not  detached  or 
reattaching  prior  to  reaching  the  tlap;  an  assumption  in 
question  when  t  low  is  analyzed  in  photos  at  these  Mach 
numbers.  Third,  that  the  '  <  w  condition  at  the  point  and 

times  in  question  are  ail  in  one  .stale,  turbulent  t  low.  Th*. 

tlow  may  transition  back  and  '  nth  depending  on  t  In.-  entiy 

conditions  such  as  angle  ot  a  t  t  a  c  k ,  d  e  t  I n >  i  i  . .  i :  .  >  n  ■  i  1  •  , 

Reynold's  number,  All  these  tact  >rs  complicated  the  .in..', 
while  si  mu!  t  aneously  it  is  '»*•  t  ug  at.t  *-mpte«l  to  use  ,  i  '• 

model  in  MKATKST  to  de*  »-rm  l  n«:  'he  resu  t  ant  tie  it  in 


In  using  the  Stanton  Number-Kntha  lpy  analysis  in  eliininat  mi>j 
dependence  on  Reynold's  number  it  produces  the  equation 


0  =  St/„VeHs  =  St(y^Vaa/2)  (  (  ** /<*)*<  )  (2  1 

where  /t,  and  Va  are  the  downstream  denr-M  '  and  velocity  and  /to 
and  Voo  are  the  fcreestrearn  density  and  velocity.  Assuming  that 
the  Stanton  number  changes  are  i  ns  iyn  i  t  i  cant. ,  this  produces 

(»'jV1//.V(0)  »  t  (•*►«(«»■ )  (22) 

and  therefore 

0  -  St^0(V«3/2  )tl«(J,r)  <2  0 

or  rearranging 

0/St  tS,(V«,3/2)  =  t(oc,  Jir)  -  O/Oret  •,/•!  • 


where 

the  new  gret: 

could  be 

<jret  = 

:  St  /^(M^/2) 

-  ( Ct/2  )  /£»  ( Voo 

3/2) 

<  2 

Using 

Spalding  and 

Chi  conditions,  St  and  Ct 

are 

t  i  Xed  a  t 

g  l  v  e  n 

reference 

value  (a.o ) . 

The  0  to  put 

into 

the  IP.  AM-. 

program  could  now  be  this  new  0  generated  by  t  h  i  s  met  Imc  vi 

0  -  t  (ch,  2>  e>  )  -  To  ♦  to»(wi  -t»<i )  +  f  i (  i  «,  -  Sllt  ’  '  'i  • 

or  0  =  Oo  +  0^<a»  -  <<  >  4  (  <(•#  -  J»f.)  <  .■  /  • 

our  basic  HKA'^KST  equation.  The  difficulty  is  we  realty 
know  the  -  t  flow  conditions  m  hir  to  the  t  1  <i 1 »  •.  i .  w  % 


using  the  treestream  conditions  which  may  not  Me  the  some. 
However,  it  was  demonstrated  earlier  at  least  tor  the  region 
ot  interest  these  conditions  cnmiii  neu  with  this  method  produce 
good  results. 


where  the  program  uses  Oret  as  the  stagnation  condition  and 
ket(l)  is  zero.  It  was  assumed  the  Q  ^  derivative  was 
insignificant.  At  the  point  in  the  entry  where  the  flight 
test  maneuvers  took  place,  HfcATKST  was  repeatedly  run  in  m 
attempt  to  tind  Q»*  and  0j4,  by  taking  that  was  assumed  t  .  ■ 

correspond  to  that  and  keynoJd's  number.  When  the  0-  and  <gs, 
derivatives  and  Qo  magnitude  could  not  converge  on  a  value  t.oi 
any  length  ot  time,  the  Q  (  ^  derivative  was  attempted,  to 

found  over  a  long  time  period  only  to  discover  the  slope  ot 
this  derivative  was  not.  linear.  So  the  approach  was  taken  to 
try  to  eliminate  or  minimize  the  0  effect..  This  is  one  oi 
the  three  derivatives  varying  with  time,  the  others  being  <  i,, 
and  Om,*  Hy  eliminating  these  two  by  subt  racting,  the  ot  hot  cm 
be  determined  since  it  is  assumed  to  he  only  a  function  ot  .» , 

J„,  and  Reynold's  number.  The  conditions  required  to  . .<  < 

relate  to  the  magnitude,  sloue,  a  ltd  correct  tunc'  i  ■  >  >  \  ■ 
variation  (  1  i  n  e  a  r  e  I  I  e  i  •  *  s  '  <  >  t  '  h<-  t  1  i  •  i  ti  t  data  result-,.  \  \  . . 


this  Reynold's  number  effect  is  established,  models  m  i  he 
H  KATE  ST  can  be  modified  and  then  other  heating  derivative.*, 
estimated  during  the  maneuvers.  It  is  desired  then  to  p'a. 
the  0 derivative  back  into  the  program  as  a  functional 
variation  that  would  make  it  linear,  a  heating  reference 
variable  such  as  0  Kckert  that  could  make  it  linear,  or  a 
reference  value  that  would  completely  eliminate  the  Reynold's 
number  derivative  such  as  in  equation  (25). 

By  assuming  we  know  Do,  0«  »  anfl  0  A  e,  ,  then  the 

g 

derivatives  can  be  found  by  solving  (when  using  Re’  ) 

He  ‘  *  0  „e  (Re  ',-Ke0'ir  )  =  O-Qo-Q^  («  -  «. ) 

"0  i  af  (  i  sr  ~  J er0  )  (  2  1 

The  program  can  evaluate  0  using  any  reference  we  wish.  when 
this  is  plotted  verses  a  Re  number  that  produces  a  'inner 
slope  the  result  should  be  a  nearly  straight  line  w  i  I  h 
bias  or  magnitude.  A  AQ  which  equals 

AO  *  D-0  (<w-«0)-o^  (  &9f-  Sgfg  )  -  (Qo+Ke  term)  (  nil 

it  now  placed  in  the  HKATKST  program  Lor  numerous  sequent  la1 

time  segments  of  approx  3  mate ly  t  i  ve  seconds  duration  each  and, 

v 

assuming  that,  the  ke  term  was  zero,  results  in  an  apw  o  : 
sloping  line  representing  the  0‘>  magnitude  and  the  Re  t  .  - 1  a 

a  specific  reference  cw0  and  Set..  Thus  the  cv  and  -‘o.  h,iv<  ').•< 
removed  leaving  the  Re  eft. -cl  and  some  tgu  magnitude.  |,,n 
computer  always  obtains  the  ■_><>  magnitude  and  :  uni  a  probit  m 

except  in  very  extreme-  c .  •  . .  •  .v  h  >• , ,  i  new  ■■  i  <  •  '  t  ■  .  •  a  '  d  '  <  ■ 
considered.  '  t  the  r  i . ,  t»  t  ?  .  ■  .  •  .  .  ■  n.*a  >  no  m-  «  i-  1  •  ■  >  ’  ■  *  •  •  -a 


be  determined  the  result  would  be  a  straight  line  or  zet o 
slope  parallel  with  the  Re  number  axis  1  ettect  eliminated) 
with  a  constant  magnitude  of  <jo.  It  the  magnitude  Qo  were 
also  subtracted  out,  the  line  would  be  at  zero.  This  result 
is  our  goal,  accomplished  by  changing  the  X  axis  Re  number 
value.  It  the  slope  is  not  zero  but  reduces  the  slope  ot  t  he 
Reynold's  number  effect,  then  it  is  still  helpful  in  out 
analysis  because  it  can  then  be  assumed  that  0  i  ■■■ 
essentially  constant  over  the  short  time  segment  in  the  eut i / 
where  a  maneuver  is  executed.  All  of  those  results  ussuiim 
that  you  have  the  correct  and  Qstr  to  subtract  out.  0^  aim 
Qr  may  be  input  as  functions  of  Re  when  maneuvers  <i  r  e 
available  at  various  Mach  numbers.  It  the  "perfect."  model  i. s 
found  the  Q  may  be  eliminated. 

From  the  previous  analysis  of  the  models,  the  ones 
showing  the  most  promise  to  be  used  in  this  approach  were  i In- 
Stanton  Number-Enthalpy  results  using  Spalding  and  chi  and  the 
Eckert  flat  plate  model  tor  turbulent  tlow.  Howevei  ,  Uu 
technique  was  first  used  to  observe  it  the  0/Uret  stamiat  ion 
model  could  be  used.  As  a  first  attempt  at  analyzing  1  he 
derivatives  and  there  effects  on  the  slope  ot  this  he.it  i>,  , 
model,  the  RC  circuit  etna  logy  was  used  tor  plotting  • 
stagnation  verses  Re'  .  'I'he  result  ot  "turning  <>•  "  tm 
derivatives  it  was  hoped  wmiln  result  in  each  re  m>  ■  v  i  !-■  '  1 " 
ft  e  a  t  i  ng  effects  they  c  re  n  I  e  ■«  no  remove  the  slopes  < .  •  >  1 1  i 

derivative  in  the  model.  if  all  t  he  derivatives  i  imm  *  n  i  . . 

the  model  were  correct  ,  the  i  '.--nil  I.  using  the  Rt  iiigi  . 

won  hi  produce  a  new  A g  w  h  i  •  •  '•  w  on  '  d  he  a  I  i  n e  i  <  *  i  . 


9 

Figure  16  displays  the  AQ  stagnation  Versus  ke  '  using  the  kC 
approximation.  The  result  indicates  that  the  slope  ot  the 
line  with  all  derivatives  input  'showing  heating  etlects  fcrom 
no  sources)  changes  at  least  once.  This  line  is  compared  to 
the  line  generated  from  the  In.-a*  my  inode I  reference  with  the 
derivative  0j<f  removed  to  include  body  flap  effects.  as  the 
comparison  indicates,  during  the  first  and  second  t Lap  sweep 
the  flap  effect  is  significantly  reduced  but  not  totally. 
Also  evident  is  that  the  slope  during  *  he  third  sweep  changes 
in  the  opposite  direction.  This  indicates  that  the  flap 
derivative  removed  may  be  too  high.  The  lower  mag.  i  tude  ot 

the  heating  in  the  line  produced  wi  t  h  is  due  to  t  lie  change 

.  ; 

in  the  body  flap  position  between  the  ke  ot  .55KS  and  I.vkk 
from  14.9  to  13.3.  Likewise,  it  the  0  derivative*  were 

removed  here,  the  ex  being  slightly  higher  (excluding  the 
maneuvers)  would  produce  a  shift  in  the  magnitude  positively. 

The  Kckert  model  was  studied  in  a  similar  approach  tot 
comparison.  Figure  17  compares  the  AQ  Kckert  verses 
with  no  heating  effects  and  with  the  0jar  derivative  removed. 
The  discontinuity  is  more  apparent  in  this  case.  The  bo.iy 
tiap  derivative  removes  most  ot  the  tlap  effects  and  changes 
the  slope  ot  the  function  in  the  same  manner  as  it  did  in  in 
stagnation  model. 

The  use  ot  the  !<C  analogy  produces  noisy  and  posstoly 
inaccurate  results.  To  !>e  .cite  the  model  can  or  c.on.oi 
used  in  HKATKST,  the  sto.in.it  urn  and  Kckert  models  were  g’.i  <-o 
into  t. he  program  with  and  without  the  I*  collection  !  :  '  i  . 

The  purpose  in  this  was  to  .*,«•••  '  t  he  y/yret  Kckert  t  .  ■»..  .  * 


AQ  €C-:ERT 


Ke  '  '  's  \  nr  PC  .  fox;,  nation 


If  it  was,  it  could  be  put  back  in  to  HEATEST  as  a  variable-  as 
described  earlier.  To  input  this  into  HEATEST  the  initial 

t 

conditions  were  established  at  the  start  of  the  Re  value  or 
interest.  HEATEST  was  allowed  to  run  tor  increments  or  r i ve 
seconds  each  generating  a  Q/Qret  Eckert  value  tor  each 
segment.  This  was  done  to  eliminate  noise  and  yener.it  e  a 
slope  from  the  combination  of  rive  se  oi.d  "steps".  The  plot  <>t 
this  run  and  a  similar  one  done  tor  0/0 ret  stagnation  are 
shown  together  in  Fig.  18.  The  results  indicate  again  that, 
the  Eckert  has  a  larger  magnitude  and  slope  compared  to  the 
stagnation  model.  Time  did  not  permit  the  insertion  or  the 
heating  derivatives  to  confirm  the  preliminary  conclusions 
drawn  from  the  RC  analogy. 

Time  did  not  permit  the  input  of  the  Stanton  Numh.-i- 
Enthalpy  model  with  the  Spalding  and  Chi  ct  values  into 
HEATEST.  It  it  were  used,  the  same  y/Oret  tor  its  mode  I  woo'd 
have  been  performed.  If  the  results  confirm  the  calculat  ion, 

if 

performed,  this  model  for  y/yret,  when  calculated  votscs  i<«; 
would  produce  a  curve  that  would  eliminate  effects  ami 

possibly  effects.  Thus  the  flight  data  would  match  the 
model  for  the  maneuvers  with  the  exception  of  changes  in  body 
flap.  It  they  were  held  fixed,  then  the  model  most  1  i  k  *-  I  y 
would  produce  an  upward  sloping  curve.  'I  li  1  s  curve  •  :ou1<* 
similar  in  appearance  to  'be  models  shown  in  Mo.  '  t'  «  1  *  '  ■ 
calculations  indicate  it  wonlu  be  much  more  linear.  '  1  ■  ■  -  o  •. , 
would  essentially  be  for  a  constant 

A  1*1  —  0  1  ‘  ^  ^  v  ' 
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B  IS 


From  AO#  a  i#F  value  can  be  input  to  produce  the  tunct inn 

0  «  O/Oretf^.^  =  0*  <x  ,  )  -Q*,,  (  iBf  -  i  )  <*.<) 

i  c  i  . 

From  the  above  equation  the  derivatives  0  bf  *  0.x  <imo  i  he 

magnitude  Qo  would  all  be  established  and,  it  done  cot  re.  '  ' y, 
should  produce  a  "match"  with  the  flight  data.  This  then 
would  be  the  new  process  ot  determining  the  heating  rate 
derivatives  tor  the  body  tlap  trom  HEATEST. 


An  attempt  to  use  the  HEATEST  procjram  in  its  oriijm.d 
torm  tailed  to  properly  determine  all  the  heating  derivatives 
that  ettect  the  body  tlap.  The  angle  oi  attack  and  body  1  l  .ip 
detlection  effects  were  approximated  but  definitive  values 
could  not  be  veritied  due  to  the  problems  in  calcul.it  ing  the 
Reynold's  number  effect.  The  Reynold's  number  effect  w.is 
found  to  be  nonlinear.  A  scientific  approach  was  used  to 
determine  it  a  model  could  be  found  that  would  linearize  <>t 
eliminate  the  Re  effect.  Several  techniques  were  studied  and 
at  least  one  was  found  that  could  be  applied  to  hkatkst.  i  iu- 
lndications  are  that  the  point  studied  on  the  flap  was  in 
turbulent  flow  but  the  possiblity  still  exists  that  it  caitd 
have  been  in  transition  or  a  local  flow  reattachment. 

The  question  still  exists  as  to  how  this  approach  ■  ,n 
be  used  on  the  centerline  body  flap  thermocouple  pond  wui<h 
appears  to  be  affected  much  more  by  elevon  detlection  than  t he 
one  studied.  The  analysis  done  in  this  report  is  only  the 
tirst  step  in  studying  the  overall  heating  effects  in  t  1 1  .  : . 
area.  It  is  a  critical  point  in  later  missions  atm  it-, 
determination  is  vital  prior  to  the  high  cross-range  m • ss i  -  ms 


in  the  future. 


Recommendations 


The  Spalding  and  Chi  technique  incorporates  many 
assumptions  but  if  used  at  a  reference  and  condition  could 
be  used  as  a  valid  Qret  mode!  in  HKATKST.  I  recommend  that  an 
attempt  be  made  to  incoporate  the  b'n  t  ha  1  py-S  t  a  nton  number 
relationship  using  the  Spalding  and  Chi  Ct  values.  This 
appears  to  be  the  most  accurate  means  ot  determiuino  t hr 
heating  effects  during  the  maneuvers. 

If  the  Spalding  and  Chi  technique  cannot  be  used,  i 
recommend  attempting  to  use  the  Kekert  turbulent  t low  eonai  ion 
or  attempt  to  use  models  where  Reynold's  number  is  not  usee. 

The  HE ATE ST  program  presently  uses  the  Kalman  f i 1 t ei  to 
take  the  equation 

X(  t )  =  F x  (  t  )  +  C.U(  t )  '  M  > 

and  produces  the  new  equation  after  the  time  increment 

X(t+At)  =  Fx(t)  +  (Jl>(t)  It.’' 

The  function  Fx(t)  has  the  value 

Kx  ( t )  *  e  At  MSI 

This  function  could  be  left  alone  after  the  first  itft.it  ion 
and  save  computer  time.  However,  in  a  related  thfs> 
student  has  taken  advantage  ot  the  A  matrix  be  i  nu  tio  c. 
and  has  subtan*  ually  r  #•  >  ’n  i'*-i  1  *  h«-  computer  time. 

I  recommend  that  as  a  means  o*  determining  t  In-  >■ .  ,  •  1  , 
n  um  be  r  that  UK  AT  K  S'!  donh:  i  >,  co  r  j >r  a  »  ••  a  1  1  aval  1  •  ’  •  1  .  '  i 
lest  Maneuver  data  . .  ' ,  and  titer  mis:.)  m:  w  1 


specific  time  interval  where  the  same  conditions  occur  in  tin 
entry  and  the  derivatives  compared  to  obtain  .1  more  acoura*. 
value  ot  . 

I  recommend  a  possible  change  in  the  flight  t*-si 
maneuver  philosophy  by  going  to  shorter  maneuvers  wh.-r. 
Reynold's  number  is  not  a  problem.  Another  approach  may  be  t. 
do  the  opposite  by  trying  a  series  ot  maneuvers  sequenced  out 
to  get  slopes  and  derivatives  ot  other  parameters  leaving  t  ti«. 
Reynold's  number  ettect  to  be  established  last. 

Due  to  the  tact  that  there  was  so  much  difficulty  1 1 
modeling  the  entry  heating  at  the  body  tiap  it  is  apparent 
that  one  or  several  things  are  takiny  place  that.  1  am  n  .< 
aware  of.  Rased  on  my  analysis,  I  would  recommend  iht' 
priority  emphasis  be  placed  in  determining  what  is  hapueniu, 
here.  More  flight  test  maneuvers  should  be  performed 
analysis.  There  should  be  maintained  a  means  of  obtaining 
thermocouple  data  from  these  and  othei  body  flap  locations  on 
future  flights.  A  mission  with  a  landing  at  Kennedy  space 
Center  or  using  properly  located  ARIA  aircraft  would  allow 
thermocouple  data  to  be  sent  in  real-time  during  the  entry 
(assuming  little  it  no  blackout  'ate.*,  placet  removing  tin*  need 
for  the  recording  and  post  Might  playback  on  t  tie  orbitet 
(  w  h  t  c  It  has  been  unsuccessful  in  se  ve  r  a  •  previous  .<  t  >  t  <n  •  • 
Until  t  tie  heating  can  tie  t  e  I  *  a  t>  1  y  predicted  in  this  1  ,  ,  • 

a  Iso  recommend  that  no  WTR  missions  he  performed  w '  * ,  •  ,, 

angle  of  attack  and  maximum  aft  CCi  condition  , 

producing  possible  excessive  ht.-at  ing. 
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Dfirljya-tiQa-.ot -H&atJLau- Model  soluti u ns 

To  perform  the  to. flowing  calculations  re  tore  nee  uoini 
were  chosen  at  various  times  during  the  point;  in  the  entry 
when  flight  test  maneuvers  were  being  performed  on  STS -J.  In 
order  to  do  the '  ca  leu  1  a  t  i  ons  the  body  flap  had  to  be  assumed 
fixed  during  these  times  and  a  second  set  of  ca  I  cm  ’  .<  t  i « <•.  s 
performed  at  one  time  in  the  entry  with  the  flap  beimi 
deflected.  The  table  of  these  reference  points  is  shown  in 
Table  2. 

Procedure  tor  Calculations  Using  Reynold's  Analogy 
Using  Reynold's  Analogy 

N 


Ct/2  «*  A/Re 

i  >  t 

St  ■  C f/2  =  Nu/RePr 

(  J  ) 

0  3  h(Taw-Tw) 

(  n 

Nu  **  hL/K 

(  A  ) 

Substituting  (3)  into  (4) 

Nu  =  QT./K  ( '’’fiw-Tw  '  '  '>  ) 

Substituting  (2)  into  (1) 

N 

Nu/RePr  -  A/R**  ■  1  1 

Substituting  (  b  )  into  <  *>  '  an  ’  r  ••  1 1  t  ang  i  in  i 


l?le  2.  Entry  Re .t e r e nce_CQn d i Li ons. 
bor  Heating  We tji.o.j  raJ.Cula L-Lons 


Ref 

Voo 

Mach 

VL.  K-7 

Loy 

Re 

T 

oc. 

n  r «;  t 

Rt 

t  t/sec 

1 

slugs/ ft 

R*» 

K6 

v1t-<  l  R 

rliuj  fc 

1 

24178 

26.1 

.59 

5.75 

.56 

35b  .7 

39.11 

6  0  .  H 

2 

22605 

24.2 

1.15 

6.00 

1.00 

36-1  .2 

39.3  l 

67.4 

3 

20985 

21.8 

1.74 

b  .  1  25 

1.33 

386.9 

38.08 

66.0 

4 

19977 

20'.  1 

2.29 

6.20 

1.59 

412.2 

44.6b 

6  1.6 

5 

18804 

18.6 

2.80 

6.25 

1.78 

425.8 

41.51 

57.6 

b 

17203 

1  6.5 

4.50 

6.39 

2.46 

446.5 

40.69 

54  .  t 

7 

15573 

14.8 

6.50 

6.50 

3.16 

462.6 

41.33 

46.8 

dt 


QL/K ( Taw-Tw )  =  APrRe 


U) 

The  result  is 

1— N 

Q/Qref  =  ( AK(Taw-Tw)/L  Re  )/Qref  < H ) 

RTo  =  (Taw-Tw)  (')) 

1  -N 

so  •  0/0 re f  =  ( AKR/LOref  )ToRe  Mm 

2 

With  To  proportional  to 

2  2 

T*/To  =  ( 1-  .  b (  -1)^)  =  (  1+  .  2^ )  Ml) 

2 

To  (  1  +  .  2M#a )  T  (  1 

with  K*.0S4  BTU/Hr-Ft  which  is  changing  and  is  higher  with 
increased  Mach  number  and  viscosity.  rile  recovery  tactor  ( R ) 
is  equal  to  .9  and  L“]07.b  feet.  In  general  K  uses 
Sutherland's  Law  of  Viscosity  which  is 

>  •’>  * 

/*  »  2.27K-8  ■{'„  /(•!•  +19U.6)  (13) 

K  =  Cp^/Pr  wil.h  Pr=,74  and  Cp-tiOUfe  (N) 

or  /<  *  ^.V^L/Re^  (  1  i 

Placing  the  known  constants  into  the  equation  the  0/tjret 

2  1  -N 

0/Ore t  =  67.9S/<AM  +  .2M^>  (T^Ke  m  > 

where  tor  laminar  flow  A-.4M  ami  N-.S.  '•'or  turbu'cni  i  ow 

A=».U296  and  N-.7.  I'simj  the  free  stream  values  at  the  ..or..':, 
of  interest,  »  he  0/0 ret  values  tor  laminar  ami  t u rbu  '  <ui  • 


were  ca  1  cm  1  .it  od. 


The  r  e  su  1  '  ,  -ire  )  ti  r.i ' .  '  •  ■ 


Table  3.  Method  I.  Calculate  ons  Using  Reynold's  An  a  J.o<j  y  and 

2.  y.t  he  rl  a  n  d  '  s_  Law 

Log  Re  Re  yU  Ref.Pt  Qret  Q/Qret  O/nr  ft 

(slugs/  Laminar  TurouJe 

fct  sec) 


5.75 

39811 

2. 74 K-7 

1 

60.8 

.02  35 

.  1  !  9  7 

6.00 

63096 

2.79K-7 

2 

67.4 

.0270 

.15]  9 

6.125 

79265 

2.95K-7 

i 

65.0 

.0294 

.  \  8  0S 

6.2 

91397 

3.11 K-7 

4 

t.  .i  #  • 

.0316 

.  20  40 

6.25 

100000 

3 . 1 8K-7 

5 

57.6 

.0333 

.2228 

6.39 

129643 

3.32K-7 

6 

54.  L 

.03  60 

.  2655 

6.50 

158489 

3. 4 4 K-7 

7 

46.8 

.04  10 

.  3  2b0 

Method  2  using  the  normal  shock  tables  required  use  or 
the  shock  tables  in  Ref.  13.  Using  the  conditions  known  ,1 t 
each  point  upstream  of  the  shock,  the  tables  allowed  the 
values  M  and  the  ratios  Ta  /T ,  and  P^/Po  to  be  found.  Using  the 
formula  from  Ret.  7  for  stagnation  pressure 

Pns  ■  P^l 1+V a/ 1 7 1 hT^ )  in  atmospheres  (17) 

we  use  this  value  in  Ref.  7  to  calulate  Log(a/uo  and  yU.  Using 
(15),  Re  is  found  and  used  with  (16)  to  find  Q/Qret  turbuien*. 
The  downstream  temperature  was  also  calculated  using  *  he 
treesteam  condition  and  the  shock  table  ratio.  The  list  « » ♦ 
results  is  on  Table  4. 

Method  3  considered  the  downstream  flow  prior  ho  the 
boat  tail  expansion  to  be  approximately  Mach  3.  Hy  annum  i  ng 
that  the  body  flap  is  tixed  at  +7  degrees  the  tables  in  *.v  '  < 

can  be  used  to  determine  a  va  lue  0 .  Using  that  d  ,  t he  imr  no  1 
Mach  number  (Mn)  can  be  found  using 

Mn  *  3  SIN  &  ( i « ) 

The  weak/strong  shock  charts  in  Ret  13  and  the  troestream 
conditions,  the  ratio  T^/To  and  thus  To  can  be  round.  With 
the  assumed  Mach  3  condition  the  rat  10  T,/To  and  r,  can 
determined  the  same  way.  The  fat  i.o  T^/t,  am:  then  rg  <■  n,  ,..w 
be  determined  using  the  previously  determined  value?,  .in  i 
tables  tor  the  Mach  number  Mn.  The  value  T*  is  now  tonne 

T*  -  T.  +  .  5  (  Tw-T  .  )  ♦  .  2  2  <  TR-T  ,  where  TK- . 7  IV,  '  • 


Table  4.  Method  2.  Calculations-  Using  -Normal-  StiocX  Tables, 
Rock  we  1 1  Normal,  stock-  Approach. 


Ref  Pt 

Mach  # 

Ta/T, 

Table 

Pns 
( atms ) 

M 

E-6 

T, 

deg  R 

Ke  ' 

O/'jrn  r 

T'l  t  r  >u  1  •  n  t. 

1 

26.1 

133.4 

.0164 

3.3 

47584 

87388 

.(>‘•7 

2 

24  .2 

114.8 

.0279 

3.2 

41811 

144175 

.  '  2  1 

3 

21.8 

93.4 

.0367 

3.17 

36117 

1 90306 

.  1  4 

4 

20.1 

87.5 

.04  34 

3.09 

36068 

2  3  3013 

.  !  68 

5 

18.6 

68.2 

.0469 

3.01 

29043 

266572 

.178 

6 

16.5 

53.9 

.0631 

2.86 

24057 

377173 

.207 

7 

14  .8 

43.5 

.0743 

2.77 

20142 

47964  L 

.247 

using  (13)  to  tind  the  viscosity  and  (15)  tor  finding  No  'he 
O/Qref  values  can  now  be  determined  us.,  ng  equation  (16).  The 
values  ot  Q/Oret  can  be  determined  the  same  way  tor  various 
changes  in  body  flap  deflection  by  fixing  all  these 
deflections  at  one  condition  (reference  point  5  was  chosen). 
This  angle  is  input  into  the  tables  and  a  new  value  ot  0  . 
The  exception  to  this  is  when  the  jangle  is  0.  Then  there 
is  no  deflection  or  shock  after  the  Mach  3  flow.  Therefore  mi» 
is  equal  to  3.  The  results  ot  these  calculations  are  listed 
in  Table  5. 

Method  4  used  a  pressure  ratio  relationship  and  • be 
angle  the  flap  makes  with  the  treestream  flow  to  determine  the 
Newtonian  flow  initial  conditions.  This  pressure  ratio  is 

2 

*VPp.  “(SIN  ^^V«Vp«)+l  (2H) 

where  &  *  (cx+J*, -7)  (  j  >  ) 

The  flap  angle  was  again  fixed  at  +7  degrees  which,  due  to  the 
boattails  slope  approximately  makes  the  flap  deflection  equn 1 
to  the  angle  of  attack.  The  shock  tables  give  us  '''he 

value  Pa  can  then  be  determined.  Using  this  value  and  the 
value  from  the  previous  method,  the  equation 

A  »  P,  /RTa 

produces  the  downstream  density.  The  downstream  Mach  num‘..-i 
can  be  found  from  the  tables  anti  velocity  Vj  can  be  exg>  i 

V,  *  Ma  a  “  (J'R'1''  ■  '  »  ' 

Reynold's  number  is  determined  using  (!6i  using  then  »  »  •  ■  ■ 


Method  3.  Weak  Shopk,  Qa  I  r».i  1 .?  t  ions  Usinj  r  i 


Table  5. 


For 

Body  Flap  Deflection  +7  deg, 

0  =  25,  Mn=l  .27  and  Ta 

/T,  =1.17 

Ref 

Ft 

T^/To 

To 

T, 

T* 

T* 

Re  ? 

6-6 

O/Qret 

1 

.0079 

45152 

16  1  24 

18847 

1  5290 

2.77  6254 

.  1  88 

2 

.0085 

43024 

15364 

18006 

14  6  80 

2.71  10238 

.  2  39 

3 

.0122 

37202 

13285 

15570 

12947 

2.54  14165 

.278 

4 

.0165 

33787 

12065 

14140 

11968 

2.44  17503 

.  O.'b 

5 

.0143 

29886 

10671 

12506 

10623 

2.30  20553 

.33! 

6 

.0180 

24806 

8858 

10382 

9566 

2.17  29256 

.  3  92 

7 

.0223 

20721 

7400 

8672 

7868 

1.96  39405 

.462 

For 

>U“18.6 

,  To=2988b ,  T, 

=  10671  , 

Tpa/To*  .01425 

Ref 

Pt 

Of  i 

Mn 

T»/T, 

** 

T* 

.  t 
Ke 

E-6 

0/0  l«‘t. 

8 

0*0  — 

— 

.357 

10670 

10109 

2.24  21004 

.  3  30 

9 

3.5  22 

1.124 

1  .078 

1  1503 

10343 

2.27  2080b 

.33' 

10 

7.0  25 

1.270 

1.172 

1250b 

10623 

2.30  20553 

.33' 

11 

10.5  28 

1  .410 

1.261 

134  56 

10890 

2.33  20366 

.  <32 

12  14.0  32  1.567  1.365  14566  11206  2.36  2013U  .'Ml 


previous  method.  Re'  and  Q/Ore£  using  (16)  completes  the 
analysis.  The  changing  body  flap  calculations  are  done  in  a 
similar  manner  resulting  in  a  changing  M4  and  other  downstream 
values.  The  resulting  Q/Qref  showed  very  little  change  with 
the  body  flap  deflections.  The  results  are  shown  in  Table  n. 
Method  5  incorporated  the  fcekert  flat  plate  formula 

Q/Qret  »  (  .0296(  RTo-Tw)  (  Re/ *)  (T*/l*  )' (24) 

*  V^Cp/778 . 3Pr  ‘ ')(P/0 

tor  turbulent  flow  with  deflection  where  T*  is  using  (l‘i)  an.t 
Re  from  (15).  All  values  used  were  the  original  f rees t ream 
conditions  at  the  reference  points,  Pr=.72,  and  Cp=6006.  ''he 
q  dot  value  is  then  divided  by  the  qret  tor  that  point.  'the 
results  are  in  Table  7. 

Method  6  used  the  previous  Newtonian  calculations  >»>■ 
and  V  with  the  enthalpy  equal  to 

Hs  -  h,m  *  V00J/2  *  CpTo  (25) 

Both  methods  to  obtain  Hs  were  used.  Then  using  a  St-.OUi  ( j 
was  calculated  from  the  equation  derived  in  Ref  13 

■  t  10  i 

0  -  ARe  ( /t  Hs/x)  -  ARe  ( JA  H  s  /  x  )  /  R  e 

*  A(/Vx^Z/ )  (  Hsxz/x  1  /  Ke  ”* 

»  .003  Va  Hs/778 . 3  •  .'*•  1 

This  was  then  divided  by  the  reference  Orel  to  obtain  '  he 
Q/Qref  values.  These  results  are  displayed  in  Table 


Table  6.  Method  4.  Newt o nian  Flow  Ana  ly sis 


For  a 

body 

flap  deflection  of  +7 

deg  and 

M  =2 .6 

Ref 

Ft 

pa/p« 

Pc 

Pa 

E-6 

a 

a 

V* 

Re  ‘ 

(.i/nro  i 

1 

380.51 

.0363 

13.82 

.43 

6738.4 

17520 

23415 

.64 

2 

324.93 

.0722 

23.45 

.76 

6578.0 

17103 

37.173 

.85 

3 

252.52 

.1158 

29.26 

1.10 

6116.5 

15903 

49428 

.06 

4 

278 .02 

.1624 

45.15 

l  .87 

5828  .9 

15155 

7506H 

1  .  3" 

5 

213.27 

.2049 

43.70 

2.04 

5481.8 

14253 

80505 

1.28 

6 

166.41 

.3450 

57  .41 

3.22 

4980.0 

12948 

11278 

1  .  *!  M 

7 

134.93 

.5133 

69.26 

4.65 

4564.8 

11868 

15436 

1  .77 

For 

changing  body 

flap# 

^=18 .6 

and 

=  .2049 

Ret 

Ft 

a  p, 

deg 

/Ppo 

Pa 

E- 

1 

6 

a  t  V  3 

■/ A 

E-6 

Re'  ^ 

0/0ret 

8 

34.5 

156.01 

31.98 

1.75 

506  3 

15190 

2.24 

76545 

1 

.20 

9 

38  .0 

184.15 

37.73 

1.91 

6297 

l  4  9  8  3 

2.27 

8226  j 

1 

.  29 

10 

41.5 

213.27 

43.70 

2.04 

5482 

14263 

2  .  10 

80505 

1 

.28 

11 

44  .0 

242.59 

49.71 

2.15 

5686 

14  0  4  9 

2 . 33 

8  2  2  6  1 

t 

.  1 4 

12 

48.5 

272.04 

55.74 

2.23 

5916 

13607 

2.36 

8  172  1 

! 

.  3 1) 

cl 


Ret  Pt. 


10.41 


15.35 


15.50 


15.39 


13.94 


14.0b 


12.41 


Q/Qret 


1726 


.  22  77 


.2384 


.2428 


.2421 


.259a 


.265  1 


Kef  Pt.  Hs-V^,  /2  0 

E8 


2.923 

8.409 

2.555 

12.767 

2.202 

14.847 

1.995 

21.800 

1.768 

19.81 

1  .480 

23.78 

1.212 

25.79 

Mil 


O/Qret 

Hs=CpTo 

E8 

0 

O/On-t 

383 

2.712 

7.802 

.  i  2  h 

.18  94 

2.584 

12.913 

.2280 

2.234 

15.060 

.2  111 

.34  30 

2.0  29 

22 . 1 67 

.  34  9l 

.34  4 

1.795 

20  .  !  1 

.349 

.440 

1.490 

23.94 

.440 

.5  50  • 

1  .245 

26.47 

.566 

Method  7  used  the  Spalding  and  Chi  technique  d  i  s<-r  i  t  j  -  -  ti 
in  Ret  11.  This  method  assumed  that  the  edge  conditions  next 
to  the  wall  and  temperature  at  the  point  ot  study  are  known. 
From  Method  3,  T2  and  Re2  are  used  as  the  edge  conditions. 
With  the  wall  temperature  known  the  ratio  Tw/Ta  is  obtained. 
From  the  tables  in  Ret  11  the  values  Fc  and  F^c  are  found. 
From  Fj  ,  the  tables  are  used  again  to  tind  Fe''t.  The  resn't 
produces  the  value  Cf  trom 

C  £  =  FcCt/FC  ( 7  ) 

Table  9  shows  the  calculated  Ct  values  and  there  input  into 
the  Stanton  Number-Enthalpy  method. 


[♦Til 


dilations  ot  Spa  'cling  and  fjhi  ivi;lm  i 
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Met.h.k, 


Tw/Te 

Rea 

K5 

Kc 

K*« 

K«*  K* 

K6 

FcCt 

Ct 

.0029 

.0059 

.0029 

.0057 

1568 

1.55 

.52 

24 

3.71 

.00  29 

.  0  0  5  6 

.0029 

.0053 

.0029 

.00  50 

26 

3.83 

.613 

1U 

3.86 

.0029 

.0047 

0/0 ret 

.  non 

.  >  HOU 
.  .1  !  • 
.ilUl) 
.  S  100 
.  M  ‘It, 


Appendix  L 


RC  Circa i t  Anal Q'JY 

An  RC  circuit,  excited  by  a  piecewise  constant  input 

yVVYV\_ 

T  j  v(t) 

\ 

It  a  voltage  e(t)  is  piecewise  constant  in  every  interval 

nT  <  t  <  (n+1 )T  (  /H  ) 

The  differential  equation  is 

RC  dV/dt  +  V  *  en  f  J'-*  * 

The  solution,  which  equals  V(nT)  at.  t  =  nT  is 

-( t-nt )/RC 

v  ( t )  *  e  n  +  [  V  (  nT  )  -  6  n  ]  e  (  n  i  ) 

Evaluating  the  solution  at  t-(n+l)T,  the  differential  e.iuut  '*m 
becomes 

-(  t-rit  1/KC 

V(nT+T)  =  en  +  (V(nT-en)le  1,1 

-T/RC  -T/RC 

V(n+1)  *  e  V(n't(  |-e  )en  ’  > 

which  is  first  order  linear.  i'dl  Mmj  en-'IVq,  where  Tea  '  >•• 
equilibrium  t  einpe »  a  t  u  re  related  t>>  *  tie  heaMnq  ininil 


e  (t) 

i 

i 

r 


Teq  =  ( Q/CCT ) 


(  33) 


and  V^T,  the  output  or  thermocouple  temperature.  The  output 
becomes 


n  -At/RC  -At/ Up  n-! 

where  T  =  (1-e  )Teq  +  e  T 


(  34  ) 


this  is  constant  at  each  At.  This  equation  could  be  used  (<> 

simulate  thermocouple  temperature  for  a  prescribed  heat  m 

corresponding  to  a  given  Teq.  Solving  for  T, 

-At/RC  -  At/  UC  n-1 

(1-e  )Teq  *  T,  -  (l-(l*e  )T  (  i‘> ' 

-At/RC  -At/RC  n-l 

or  T,  ■  Tn/ll-e  )  —  1 1/ (1-e  )  —  1 )  T  (3t>) 


rear rang l ng 


-At/RC 


-  At/RC 


Teq(tn)  «  (  l/l-e  )T  -e  T  (tn-1) 


(  i7  1 


i  t 


-At/RC 
g  -  1/ ( 1-e  ] 


(  iM  ) 


then  Teq(t.i)  ■  g  T1  -  (g-1  )Tl  ( tn-1  )  (  H  > 

This  is  the  circuit  analogy  used  to  calculate  an  equi  librium 
temperature.  It  y-l,  then  you  have  just  the  equilibrium 
assumption.  The  maximum  value  of  g  is  assumed  to  be  l.'>  foi  I 


second  ot  time  and  the  limit  of  t  is  Rcx.l. 
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